Ba + Quadrupole Polarizabilities: Theory versus Experiment 



B. K. Sahoo * 

Theoretical Physics Division, Physical Research Laboratory, Ahmedabad-380009, India 



(N 

o 

(N 
& 

in 



43 



o 

-I— > 
O 



>: 
ov 

in 

00' 

q; 

ov 
o: 

(N 



X 



B. P. Das 

Theoretical Astrophysics Group, Indian Institute of Astrophysics, Bangalore-560034, India 
(Dated: Received date; Accepted date) 

Three different measurements have been reported for the ground state quadrupole polarizability 
in the singly ionized barium (Ba + ) which disagree with each other. Our calculation of this quantity 
using the relativistic coupled-cluster method disagrees with two of the experimental values and is 
within the error bars of the other. We discuss the issues related to the accuracy of our calculations 
and emphasize the need for further experiments to measure the quadrupole polarizability for this 
state and/or the 5D states. 

PACS numbers: 32.10. Dk, 32.70.Cs, 31.15.ap, 31.15.ag 



I. INTRODUCTION 

High precision studies of polarizabilities of atomic sys- 
tems are of interest in a number of different problems 
in physics PHI]- A wide variety of methods have been 
used to calculate the polarizabilities of these systems @- 
9] . Calculations of polarizabilities depend on the matrix 
elements between different atomic states and the excita- 
tion energies between them. Therefore, the accuracies of 
these calculations depend on the uncertainties in both the 
quantities. In general, it is very challenging to minimize 
these uncertainties. However, in a sum-over-states ap- 
proach [lol , the major uncertainties in evaluating the 
polarizabilities can be reduced by using the experimen- 
tal energies. Furthermore, the accuracies of the matrix 
elements can also be improved by matching the results 
of the lifetimes/branching ratios of atomic states from 
sophisticated many-body calculations and high precision 
measurements. For cases where the experimental results 
are not sufficiently accurate or if all the available mea- 
sured results are not in agreement, it is not possible to 
test the accuracies of the ah initio results. It might be 
useful in such situations to employ the sum-over-states 
approach to evaluate the polarizabilities. 

For the ground state quadrupole polarizability in Ba + , 
the available experimental results [r3 - [l4| are neither in 
agreement with each other nor with the calculations that 
are reported in this work. Investigation of various prop- 
erties using relativistic many-body methods for high pre- 
cision studies in this ion are useful in the context of a pro- 
posed parity nonconservation experiment [TEl [l6| , search 
for the nuclear anapole moment [TtJ, estimation of the 
uncertainties for a proposed optical clock [l8|, determi- 
nation of the nuclear octuple moment [l9[ etc. In this 
paper, we report the results of our calculations of the 
matrix elements between different atomic states of Ba + 
using the relativistic coupled-cluster (RCC) method. We 



also intend to test the accuracies of some of the important 
matrix elements by using them to estimate the lifetimes 
of the 5D states and comparing with their corresponding 
measurements. 

Before presenting our results, we define the quadrupole 
polarizability for a general atomic state in the following 
section. We give a brief description of the method of cal- 
culations of the wave functions and the matrix elements 
in Section IV and then present our calculated results and 
discussions after which we make our concluding remarks. 



II. THEORY 



The potential energy of an atom in the presence of a 
static electric field is given by (for example, see Q) 



V(r) = -Pi£i - -QijdiEj 
= eY J ^r l P l (cos6), 



(1) 



where pi and Qtj are the components of the electric dipole 
and quadrupole operators, respectively, Ei and Pi(cos9) 
are the applied electric field and Legendre polynomial, 
respectively, with component i. The quadrupole polariz- 
ability of an atomic state N^) is related to the energy 
shift by the equation [2(| [2l| 

1 



5E n (£) 



- n f 2 



(2) 



where is known as quadrupole polarizability of the 
state \^ n ) which is defined as 



= -2E 



K*„iQi* m >r 



E,, 



E n 



(3) 
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and Q= q is the total electric quadrupole operator 
whose single particle reduced matrix element is given by 

(K f \\q\\ Kl ) = { Kf \\C^\\ Ki } 



drr 2 (Pf(r)Pi(r) + Q f (r)Qi(r)), (4) 
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where P(r) and Q(r) represent the radial parts of the 
large and small components of the single particle Dirac 
orbitals, respectively. The reduced Racah coefficients are 
given by 



jf k ji 



1/2 -1/2 



■K{l KJ ,k,l Ki ), (5) 



with 



7r(Z, m, Z') 



1 for I + m + I' = even 



otherwise. 



(6) 



III. METHOD OF CALCULATIONS 

In the present study, the atomic states of Ba + that 
we have considered have closed-shell cores and a valence 
electron v. We adopt a two step procedure to calculate 
the wave functions for these states in the Fock space rep- 
resentation: First, the Dirac- Fock (DF) wave function for 
the common closed-shell core [5p 6 ] (denoted by |$o)) is 
calculated. In the next step, we append a corresponding 
valence orbital v to obtain the new DF wave function for 
the new configuration (denoted by |$„) = aJ|$o))- The 
atomic state for the new configuration is expressed using 
the RCC ansatz as 



\l + 54|$„), 



(7) 



where T and S v are the excitation operators that take 
into account correlation effects arising from the core and 
valence electrons, respectively. We consider only the 
singly and doubly excited configurations from the DF 
wave functions by approximating the T and S v opera- 
tors as 



T = Ti + T 2 and S v = Sx v + S 2v . 



(8) 



The above approximation is known as the coupled cluster 
singles and doubles (CCSD) method. Also, we construct 
triple excitation configurations from perturbatively 
in the spirit of the CCSD(T) approach 0, [TJ, HHHl- 

The amplitudes for both the T and S v operators are 
determined using the following equations 



<S>«\{H N e T }\1> Q ) = 6 , K AE C . 



(9) 



and 

{*Z\{H^ T }{1 + S V }\* V ) 



($Z\1 + S V \$ V ) 

($„|{i^ T }{l + 54|$„) 
= [6 v , K + {<i>«\S v \<S> v )]AEf\ 

(10) 

respectively. Here K = 1, 2--- represents the singly, 
doubly etc excited configurations with respect to their 



corresponding reference states, i/jve T denotes the con- 
nected terms of the normal order Dirac-Coulomb (DC) 
Hamiltonian (i/jv) with the T operators, AE corr and 
AE% tt are the correlation energy for the closed-core and 
attachment energy of the valence electron v, respectively. 

We evaluate the scalar polarizabilities by expressing 
them as the sum of three contributions 



a v 2 (c)+a v 2 (cv)+a v 2 (v), 



(11) 



where a v q (c) takes into account the contributions of the 
core orbitals, a v q {cv) and a v q (v) are the core- valence and 
valence contributions, respectively. In terms of the re- 
duced matrix elements, the expressions of these parts are 
given by 



and 



al(c) 
a v 2 (cv) 



2 v \(j P \n\\Jc)Y 

5 ^ E D -EL 



5(2J„ + 1) 



\(J V \\Q\\J C )\ 
E c — E v 



2 \{J v \\Q\\J m )\ 



E v 



(12) 
(13) 

(14) 



where (J P ||Q||J C ) are the reduced matrix elements be- 
tween the atomic states with angular momenta J p and 
J c . a 2 (c) and a 2 (cv) have been calculated using the 
third order many-body perturbation theory (MBPT(3) 
method) from an expression 



a 2 l = (* n |Q|*W>, 



(15) 



where \^>n ) is like a first order perturbed wave function 
and it is obtained by solving the following inhomogeneous 
equation 



(H-E n )\*U) = (i#>-Q)|tt n ), 



(16) 



with En^ — (^nlQIv&n); which is similar to the first or- 
der perturbation equation. It should be noticed that un- 
like for the dipole operator [t| [27]], En results are fi- 
nite for the quadrupole operator and correspond to the 
quadrupole moments of the respective states \^ n )- 

Contributions from a 2 (v) were determined by calcu- 
lating important intermediate states explicitly using the 
CCSD(T) method. The reduced matrix elements be- 
tween different states were computed using the following 
expression 



(JfWQU) = 



(JtWii + sUQii + SiWi) 



(17) 



where Q 



e Tt Qe T 



and J\f v 



(* v \e 



Sje 1 " 1 e T 'S v \<& v ) involve two non-truncating series in the 
above expression. The details of the calculations of these 
terms are discussed elsewhere || [IH, HH, [23| . 
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TABLE I: The E2 matrix elements and the ground state 
quadrupole polarizability in Ba + in au. Results are given 
up to significant digits. Possible uncertainties in the results 
are given in the parentheses. 



Transition 


Amplitude 


a E2 


a 6 2 s (v) 






6si/ 2 -¥ 5d 3/2 


12.76(5) 


1466(11) 


6^3/ 2 


16.58(12) 


263(4) 


7^3/2 


5.727(7) 


24.07(6) 




4.036(5) 


10.72(3) 


5d 5 /2 


15.99(8) 


1978(20) 


6^5/2 


19.99(20) 


380(8) 




7.024(9) 


36.1(1) 


8^5/2 


5.022(5) 


16.59(3) 


al s (c) 




46(2) 


a 6 2 "(cv) 




-0.001(0) 


afttail) 




50(10) 


af (total) 




4270(27) 


Expt [12] 




2050(100) 


Expt [13] 




2462(361) 


Expt [14] 




4420(250) 


Theo [32] 




4091.5 


Theo [33] 




4182(34) 


Theo [37] 




4821 



IV. RESULTS AND DISCUSSIONS 

Below we present the quadrupole polarizabilities for 
the ground and the first two excited D states in Ba + . 
We have used the experimental energies from the Na- 
tional Institute of Science and Technology (NIST) data 
base 28j in the sum-over-states approach to evaluate the 
major contributions to the quadrupole polarizabilities, 
which are the the valence correlation effects. The main 
purpose of doing this is, as stated in the first section, 
to minimize the uncertainties in the calculated results so 
that we shall be able to compare them meaningfully with 
the available experimental results. Contributions from 
the higher excited states that cannot be accounted by the 
valence correlation in the the sum-over-states approach 
are evaluated using the MBPT(3) method; since their 
contributions are typically smaller. The uncertainties in 
these results are estimated by scaling results obtained 
using this method with the CCSD(T) calculations. 

In Table HI we present the results of our quadrupole 
polarizability calculations of the ground state of Ba + . It 
is clear from this table that the dominant contributions 
(approximately 80%) come from the 5D states. Thus the 
accuracy of the ground state quadrupole polarizability 
calculation depends primarily on the accurate determina- 
tion of the E2 matrix elements of the 6S1/2 — » 5Z? 3 / 2 and 
QS1/2 — ► 5-D 5 / 2 transitions. These matrix elements were 
also calculated by us earlier using the same CCSD(T) 



method but with different basis functions [lg, [29J , and all 
the results are in good agreement. There are also other 
calculations available for these matrix elements using dif- 
ferent variants of the RCC methods and basis functions 
[30l - [33l ] ; all the calculated results seem to be in reasonable 
agreement with each other. Moreover, the accuracies of 
these matrix elements can be verified by using them to 
estimate the lifetimes of the 51? states and comparing 
them with the measurements as stated in the Introduc- 
tion. Possible transition channels from the 5D 3 / 2 state to 
the ground state are due to the Ml and E2 multipoles. As 
we have shown in our earlier work (29| . the Ml transition 
probability is very small in this case and its contribution 
to the lifetime of the 5-D 3 / 2 state is negligible. Therefore, 
neglecting the Ml contribution and using our calculated 
E2 matrix element of the 6S1/2 — > 5_D 3 /2 transition, we 
obtain the lifetime of the 5-D 3 /o state as 79.8(6) s. The 
measured values are 79.8(4.6) s [H and 89.4(15.6) s [Hj]. 
Our result is in good agreement with the first experimen- 
tal result and is within the error bar of the second result. 
To calculate the lifetime of the 5D 5 / 2 state, it is necessary 
to take into account all the Ml and E2 transition prob- 
abilities from this state to the 6S1/2 and 5D 3 /2 states. 
We have also shown in Ref. [29| that only the E2 and 
Ml transition probabilities of the 5D 5 / 2 — > 6iS 1 / 2 and 
5-D 5 / 2 5-D 3 / 2 transitions, respectively, are significant 
in the determination the lifetime of the 5D 5 / 2 state. In 
the present work, we have calculated the above Ml tran- 
sition amplitude to be 1.544(1) au which is in agreement 
with the result reported in Ref. [29(. Using this value 
we find the lifetime of the 5£> 5 / 2 state to be 29.8(3) s 
with 84% branching ratio to the 5D5/2 —> &S\/ 2 tran- 
sition. This is also in agreement with the experimental 
results which are reported as 34.5(3.5) s [11] , 31.6(4.6) 
s [35[ and 32(2) s [36]. We have used the experimental 
wavelengths to determine the lifetimes of the 5D states 
in order to verify the accuracies of the calculated E2 ma- 
trix elements. This analysis suggests that our calculated 
E2 matrix elements are accurate to within a few percent 
and therefore they can be considered for high precision 
studies of the quadrupole polarizabilities of Ba + . 

The next significant contributions to the ground state 
quadrupole polarizability come from the 6D states. It 
is difficult to estimate the accuracies of the E2 matrix 
elements of the 6S , 1 / 2 — > 6£> 3 / 2 and 6S , 1 / 2 — > 6D 5 / 2 tran- 
sitions from the measured lifetimes of the 6D states be- 
cause of their negligible roles in the theoretical determi- 
nation of these lifetimes. Also, the contributions from the 
higher excited states (tail) and core correlations to the fi- 
nal result of the quadrupole polarizability of the ground 
state in Ba + are non negligible. A suitable method to 
test the validity of all these contributions is to compare 
the final calculated result with the available measure- 
ments. The reported experimental results for the ground 
state quadrupole polarizability of Ba + are 2050(100) [Hj], 
2462(361) [H and 4420(250) [3 in au. The first two re- 
sults agreewith each other, but they are completely in 
disagreement with the latest result. All the reported ex- 
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TABLE II: The quadrupole polarizabilitiy of the 5D 3 /2 state 
of Ba + in au. Uncertainties are given in the parentheses. 



Transition 


Amplitude 


a E2 


5d3/2/ \ 

a 2 (v) 






5d 3/2 -4 6s 1/2 


12.76(5) 


-733(6) 


7*1/2 


4.882(5) 


13.96(3) 


8S1/2 


1.558(2) 


1.0(2) 


9si/ 2 


1.360(1) 


0.668(2) 


6d 3 /2 


8.36(4) 


37.3(5) 


7d 3 /2 


3.021(3) 


3.65(1) 


8d 3 /2 


2.183(2) 


1.69(1) 


54/2 


6.83(2) 


1278(7) 


64/2 


5.349(2) 


15.21(1) 


7rf 5 /2 


1.978(2) 


1.56(1) 




1.431(2) 


0.727(3) 


5ff7/2 


8.45(8) 


26.9(5) 


6ff7/2 


8.63(9) 


25.7(5) 


5d3/2, \ 

"2 (c) 




46(2) 


5d3/2, \ 

a 2 [cv) 




-0.49(3) 


5d3/2/, 

q 2 (tail) 




116(30) 


a 5d3/2 ^ Q ^ a ^ 




835(32) 



perimental results have relatively large error bars. Two 
of these experimental results are just half of our calcu- 
lated result. In such a situation, it will not be possible to 
test the accuracies of the many-body methods that have 
been used to perform the calculations as well as those 
that are are likely to be developed in the ab initio frame- 
work to calculate these quantities. Also, the above ex- 
perimental techniques could single out the contributions 
from the 5D states to the ground state quadrupole polar- 
izability and their values have been reported as 1562(93) 
[H, 2050(100) 0, 1828(88) QJ] and 1524(8) [13 in au. 
On the other hand using the E2 matrix elements of the 



6S1/2 -> 5£> 3 /2 



and 65*1/2 ~> 5D 5 / 2 transitions, whose 



'5/2 



accuracies have been discussed above, we obtain the com- 
bined contributions of the 51? states to the ground state 
quadrupole polarizability in Ba + as 3444(23) au with in- 
dividual contributions as 1466(11) and 1978(20) in au 
from the 5-D 3 / 2 and 5D 5 / 2 states, respectively. It, there- 
fore, appears that the extracted values of these contri- 
butions of the 5D states to the ground state quadrupole 
polarizabilities in the above experimental analysis may 
not be the contributions from both the 5D states, but 
rather from either the 5D 3 / 2 or the 5D 5 / 2 state individ- 
ually in the different experiments. 

To the best of our knowledge, three calculations of the 
ground state quadrupole polarizability in Ba + have been 
carried out, yielding 4091.5 H3, 4182(34) [H and 4821 
37j in au. All these results are in agreement with our 
calculation; the first two are obtained by the linearized 
RCC method (SDpT method), in contrast to our the non- 
linear approach, but with the same level of hole-particle 
excitations from the DF wave functions. The two results 



differ from each other by about 2%, even though the same 
SDpT method was employed in the two cases, and they 
are at the lower limit of the latest experimental result 
4420(250) au. Their calculated results are slightly lower 
than our result 4270(27) au and also, they do not over- 
lap within their respective predicted uncertainties. The 
main reason for the difference between the SDpT and 
our results is that the calculated values of the E2 matrix 
elements of the 6£i/ 2 — > 5D 3 / 2 and 6Si/ 2 — > 5D 5 / 2 tran- 
sitions for the two cases are different . It is necessary 
to emphasize that the precision of the theoretical deter- 
mination of the lifetimes of the 5D states depend criti- 
cally on the accuracies of these matrix elements [2914331 ] . 
A precise measurement of the quadrupole polarizability 
of the ground state in Ba + could test the accuracies of 
these E2 matrix elements that are calculated by different 
many-body methods. A third calculation of the ground 
state quadrupole polarizability has reported a result of 
4821 au [13] using rather simple wave functions based on 
the asymptotic behavior and the binding energy of the 
valence electron in contrast to other calculations which 
are based on all order perturbative methods in the RCC 
framework. Nevertheless, all the calculated results sug- 
gest that the ground state quadrupole polarizability in 
Ba + is approximately 4200 au and all the E2 matrix el- 
ements of the 6S 1 / 2 — > 5£> 3 /2 and 6S 1 / 2 — > 5D 5 / 2 tran- 
sitions arc in reasonable agreement with each other. In 
view of the discrepancies between the calculated and ex- 
perimental results as well as among the individual values 
of the latter, it would indeed be desirable to perform pre- 
cise measurements of the quadrupole polarizability of the 
ground state of Ba + . 

An alternative approach to resolve this problem would 
be to measure the scalar quadrupole polarizabilities of 
the 5D states with high precision, and attempt to identify 
the individual contribution of the the 6S state, perhaps 
in a manner similar to the method used in obtaining the 
contributions of the the 5D states to the quadrupole po- 
larizability of the ground state. If that could be achieved, 
then it be possible to extract the E2 matrix elements of 
the 6S1/2 — ► 5D3/2 and 6Si/ 2 — > 5D 5 / 2 transitions. The 



advantages of measuring the scalar quadrupole polariz- 
abilities of the 5D states in addition to the ground state 
could be of two fold: (i) the core-correlation effects, which 
is one of the significant contributions is the same for the 
states that we have considered, can be ignored in the es- 
timation of the accuracies of the E2 matrix elements and 
(ii) By combining the quadrupole polarizabilities with 
the lifetimes of the 5D states, it would be possible to 
extract the E2 matrix elements of the above two transi- 
tions. For this purpose, we also calculate the quadrupole 
polarizabilities of the 5D states. 

In Table m we present the scalar quadrupole polariz- 
ability of the 5D3/2 state. The largest contribution to 
this quantity comes from its fine structure partner fol- 
lowed by the 6S state, but with opposite sign resulting 
in a strong cancellation between them. The other sig- 
nificant contribution is from the higher excited g states 
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TABLE III: The quadrupole polarizability of the 5D 5 / 2 state 
in Ba + (in au). Estimated uncertainties to the results are 
mentioned in the parentheses. 



Transition 



5d5/2 



5fi 



(v) 



■5/2 



5d5/2 
5d5/2 



5d5/2 



(cv) 
(tail) 



5d5/2 



(total) 



Amplitude 



6Sl/2 


15.99(0) 


7S1/2 


6.409(7) 


oSl/2 


1 QQO/0\ 


9S1/2 


1.723(2) 


54/2 


6.83(2) 


64/2 


5.754(5) 


7d 3 / 2 


2.030(2) 


8^3/2 


1.457(2) 


6d 5 /2 


11.26(2) 


7rf 5 /2 


4.029(3) 


84/2 


2.925(3) 


5#7/2 


2.96(5) 


637/2 


3.00(6) 


539/2 


10.48(8) 


6ff9/2 


10.62(9) 



-659(6) 
16.39(3) 
1.109(3) 
0.724(2) 
-852(5) 
12.03(3) 
1.114(2) 
0.509(1) 
45.8(2) 
4.38(1) 
2.05(1) 
2.24(8) 

2.1(1) 
28.0(4) 
26.3(5) 

46(2) 
-0.48(3) 
121(35) 

-1201(36) 



lying g states are significant, and it would be appropriate 
to use a method that implicitly takes all possible inter- 
mediate states into account, perhaps similar to Ref. (39j . 
for an accurate evaluation of the scalar quadrupole po- 
larizabilities of the 5D states. Nonetheless, the present 
study captures several classes of important correlation 
effects and it would be useful in guiding experiments on 
quadrupole polarizabilities in Ba + . 



V. CONCLUSION 

We have carried out a detailed analysis of the calcu- 
lated and experimental results of the reported ground 
state quadrupole polarizability in Ba + and highlighted 
the disagreement between the different studies. The re- 
ported experimental results are not reliable enough to 
test the validity of the calculated results. On the ba- 
sis of different physical considerations, we propose new 
theoretical and experimental studies of the quadrupole 
polarizabilities of the 5D states to test the accuracies of 
the E2 matrix elements between the ground state and the 
5D states and the quadrupole polarizability of the ground 
state. We have also presented the results of our calcula- 
tions for the quadrupole polarizability of the 5D states 
in Ba + using the E2 transition amplitudes obtained from 
our calculations. 



(given as tail) , which are not taken into account explicitly 
in the sum-over-states approach. Also, we find the trend 
of the contributions to the 5D 5 / 2 state scalar quadrupole 
polarizability, given in Table IIII| is similar to those for 
the 5D 3 / 2 state but, the contributions from this state and 
6S state have same sign. The contributions from the high 
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